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ABSTRACT 

We present VLT/FLAMES spectroscopic observations (R^6500) in the Ca II triplet region for 470 probable 
kinematic members of the Sculptor (Scl) dwarf spheroidal galaxy. The accurate velocities (±2 km s _1 ) and 
large area coverage of Scl allow us to measure a velocity gradient of 7.6^2 2 km s_1 deg^ 1 along the projected 
major axis of Scl, likely a signature of intrinsic rotation. We also use our kinematic data to measure the 
mass distribution within this system. By considering independently the kinematics of the two distinct stellar 
components known to be present in Scl, we are able to relieve known degeneracies, and find that the observed 
velocity dispersion profiles are best fitted by a cored dark matter halo with core radius r c = 0.5 kpc and mass 
enclosed within the last measured point M(< 1.8 kpc) = 3.4 ± 0.7 x 10 8 M Q , assuming an increasingly 
radially anisotropic velocity ellipsoid. This results in a mass-to-light ratio of 158±33 (M/L) Q inside 1.8 kpc. 
An NFW profile with concentration c = 20 and mass M(< 1.8 kpc) = 2.2iJ;° x 10 8 M is also statistically 
consistent with the observations, but it tends to yield poorer fits for the metal rich stars. 

Subject headings: Local Group — galaxies: dwarf — galaxies: Sculptor Dwarf Spheroidal — galaxies: kine- 
matics and dynamics — dark matter 



1. INTRODUCTION 

Recent wide area photometric and spectroscopic surveys 
of large numbers of stars in Local Group dwarf spheroidal 
galaxies (dSphs) have resulted in major progress in our 
knowledge of these objects and a more complete picture of 
their characteristics. They appear to be the most dark mat- 
ter (DM) dominated objects known to date, with mass-to- 
light (M/L) ratios larger than 100 (e.g.. iKleyna et al]l2005t 
IWalker et al.l 120071 and references therein); their kinematic 
status is generally believed to be undisturbed by interac- 
tions wit h the Milky Way (M W), with the exception of Ca- 
rina (see Muno zet al.J)2006L and references therein), Bootes 
(e.g., [B elokuro v et al.| 120061) . and perhaps Ursa Minor (e.g., 
iBellazzini et al.l 120021) . In some cases a number of spa- 
tial and kinematic substructures have been observed (e.g., 
IWilkinson et d]l2004t IColeman et af]l2005t iBattaidia 1 120071 
hereafter B07), which are perhaps the remnants of accreted 
subunits (stellar clusters or smaller galaxies). Furthermore, 
the combination of metallicity ([Fe/H]) and kinematic infor- 
mation have unveiled the presence of multiple s tellar compo- 
nents in the Sculptor (Scl) and Fornax dS phs dTolstoy et al.1 
12004 hereafter T04; iBattaglia et aT1l2006l) . where the metal 
rich (MR) stars have a more centrally concentrated, less ex- 
tended spatial distribution, and show colder kinematics than 
the metal poor (MP) stars. 

The presence of multiple stellar populations may require a 
modification in the way these systems are dynamically mod- 
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eled to derive their mass content. Traditionally, their mass 
distribution has been derived from the Jeans analysis of the 
line-of-sight (l.o.s.) velocity dispersion obtained considering 
all stars as a single component embedded in an extended DM 
halo. This analysis is subject to the well known degeneracy 
between the mass distribution and the orbital motions of the 
individual stars in the system (mass-anisotropy degeneracy), 
which prevents a distinction amongst different DM profiles. 

In this Letter we show that it is possible to partially break 
this degeneracy by modeling Scl as a two-(stellar)component 
system embedded in a DM halo. As part of the modeling of 
Scl, we also analyzed its dynamical status, and for the first 
time, we determine the presence of (statistically significant) 
intrinsic rotation in the stellar component of a dSph. 

2. OBSERVATIONS 

We acquired photometric data from the ESO/2.2m Wide 
Field Imager (WFI) at La Silla, and spectroscopic follow-up 
data in the Ca II triplet (CaT) region from FLAMES at the 
VLT. Our spectroscopic targets were chosen from the WFI 
imaging to have the colors and magnitudes consistent with 
Red Giant Branch (RGB) stars. We acquired data for 15 dif- 
ferent fields in Scl (7 fields in addition to those presented 
in T04). For details of the data reduction and analysis see 
Battaglia et all (120061 12008). From the imaging we obtained 
the surface density profile of the different stellar populations 
in Scl and from the spectroscopy we measured line-of-sight 
(l.o.s.) velocities, and CaT metallicities (ITe/Hl) calibrated 
using the relation derived in Battaglia et al. (2008|). 

Using an iterative procedure we measure a systemic ve- 
locity Vhoi,sys = 110.6 ± 0.5 km s _1 and a dispersion a — 
10.1 ± 0.3 km s _1 for Scl. A simple kinematic membership 
selection (stars with heliocentric l.o.s. velocity within 3<r from 
Whci,svs) produces 470 probable Scl members, as shown in 
Fig.Qji. We use this sample when investigating the kinematic 
status of Scl (Sect. [3), and adopt a more accurate membership 
criterion when deriving the l.o.s. dispersion profile (Sect. |4j>. 

3. KINEMATIC STATUS AND ROTATION 
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We first focus on the kinematic status of Scl in terms of 
tidal disruption and intrinsic rotation, as both these factors 
can inflate the observed l.o.s. velocity dispersion profile. 

No signs of tidal disruption such as tidal tails and S- 
shaped contours are found in our photometric data-set 
(Battaglia et al., in prep), although they are not nec- 
essarily expected as shown in N- body simulations by 
iMunoz. Maiewski & Johnston I d2008l) . The presence of a 
velocity gradient along the direction of the galaxy's orbital 
motion (indicated by the proper motio n direction) is an- 
other possible sign of tidal disruption (e.g. lOh. Lin & Aarseth I 
1995). We will call such a velocity gradient "apparent rota- 
tion" to distinguish it from intrinsic rotation. Fig. [TJ) shows 
the Galactic standard of rest (GSR) velocities of the prob- 
able members of Scl found within ±0.15 deg from its pro- 
jected major axis (position angle PA= 99°). At a pro- 
jected radius R ~ —0.7 deg the median velocity is 8.7i"f'5 
km s _1 larger than the systemic and at R ~ +0.5 deg it is 
4.7j^34 km s -1 lower. The best-fitting straight-line gives a 

velocity gradient of 7.6^2 2 ^ m s_1 deg -1 . For comparison, 
no velocity gradients are found along the minor axis (PA= 9°) 
and other two explored axes at inte rmediate PA (45°, 135°). 
Other techniques such as bisection dWalker et alJl2006h sup- 
port these results. This signal is unlikely to be significantly 
affected by the foreground since at a projected distance from 
Scl's center of 0.5-0.7 deg the expected fraction of MW con- 
taminants in the 3er membership range is only ~40%. 

If the measured velocity gradient is caused by tidal disrup- 
tion, then it should arise along the direction of Scl's proper 
motion, which appe ars not to be the case acco rding to the 
measu rements by ISchweitzer et all d 19951) and iPiatek et"aT] 
(2006). Although the errors on the proper motions are large, it 
is encouraging that both works agree on the amplitude along 
the North-South direction, indicating that it is unlikely that 
more accurate measurements will lead to an orbit aligned with 
the major axis of Scl. In Fig. la we also show the resulting 
orbits obtained by integrating these proper motions and the 
l.o.s. velocity of Scl in a Galactic potential. If the tidal debris 
is aligned with the orbit, then approaching velocities would 
be detected on the East side of the galaxy, and receding on 
the West side, op posite to what is observed. Furthe rmore, 
the simulations bv lMunoz. Maiewski & Johnston I (120081) also 
show that velocity gradients due to tidal disruption usually ap- 
pear at larger projected radii, equal to the King limiting radius 
(~ 1.3 deg for Scl, see Irwin & Hatzidimitriou 1995) with 
amplitudes of <2 km s _1 , much smaller than measured here. 
From these arguments we conclude that it is more likely that 
the detected velocity g radient is due to intr insic rotation than 
to tidal disruption (see lWestfall et al.ll2.Q06l for a different in- 
terpretation based mostly on photometric data). This would 
therefore be the first time that statistically significant rotation 
has been found in the stellar component of a dSph. 

In the following analysis we use rotation-subtracted GSR 
velocities. The rotational velocity subtracted to each star is 
v m t — kx, where k — —7.6+22 km s_1 deg -1 and x is the 
projected abscissa on the sky. 

4. THE MASS CONTENT OF THE SCULPTOR DSPH: TWO STELLAR 
COMPONENTS MODELING 

Scl hosts two stellar populations with distinct spatial, 
metallicity and kinematic distributions (T04 and references 
therein). Our spectroscopic sample of RGB stars shows that 
MR ([Fe/H]> —1.5) stars are less extended spatially and have 
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FIG. 1. — a) Velocity field, smoothed with a median filter, for 3cr Scl 
kinematic members. We use the velocities in the GSR frame to avoid spu- 
rious gradients introduced by the component of the Sun and Local Stan- 
dard of Rest motion's along the l.o.s. to Scl. The color bar gives the ve- 
locity scale. The dotted ellipses are placed at 0.2, 0.5, 0.7 deg to give an 
idea of th e distance scale. The solid elli pse shows the nominal tidal ra- 
dius (from Irwin and Hatzidimitriou 1995). The dashed and solid lines in- 
di cate, respectively, Scl's orbit calculat ed using the proper motion measured 
by Schweitzer et al. 1 1995) and Piatek et al. 1 2006), where red and blue shows 
the receding and approaching velocities, b) Median «gsr P er ^ m versu s dis- 
tance along the projected major axis (diamonds with error-bars) for 3a Scl 
kinematic members located within ±0.15 deg from this axis. Each bin con- 
tains ~30 stars. The solid and dashed lines show, respectively, Scl systemic 
velocity and the best fit to the observed velocity gradient (see text). Positive 
distances correspond to the east side of the galaxy, negative distances to the 
west side. 

colder kinematics than MP ([Fe/H]< —1.7) stars 7 . MR and 
MP RGB stars appear to trace, respectively, the Red Horizon- 
tal and Blue Horizontal Branch populations (RHB and BHB, 
respectively), as shown in Figure[2] The number surface den- 
sity profile of RGB stars from our photometry is well ap- 
proximated by a two-component fit, where each component 
is given by the rescaled best-fitting profile derived from our 
photometry for the RHB and BHB stars. 

4.1. Observed Velocity Dispersion Profiles 

In the following we derive the observed velocity disper- 
sion profiles of MR and MP RGB stars and use them jointly 
to constrain the Scl mass content. However, we first need 
to take into account the presence of MW interlopers which 
might contaminate our sample and affect the dynamical mod- 
eling. The probability of a MW interloper within the 3er kine- 
matic sample increases with the projected distance from the 
Scl center as a consequence of the decreasing ratio of the stel- 
lar density of Scl versus the MW. As we cannot directly weed 
out these interlopers using our photometry or spectra we take 
them into account in a statistical way using a maximum like- 
lihood method. 

The probability P(Uj) of observing a star with a velocity Vi 
and error <7j at a projected distance from the center Ri< R < 
R2, where Ri and R2 are a generic inner and outer radius 

7 We exclude the region —1.7 <[Fe/H]< —1.5 to limit kinematic con- 
tamination between the two populations. 
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FIG. 2. — Number surface density profile of RGB stars in Scl from 
ESO/WFI photometry (squares with error-bars) overlaid to the best-fitting 
two component model (solid line) given by the sum of a Sersic (dotted line) 
and Plummer (dashed line) profiles. These are obtained from the rescaled 
profiles that best fit, respectively, the distribution of RHB and BHB stars 
(diamonds and asterisks with error-bars, respectively) in Scl. The Galactic 
stellar contamination has been subtracted from each point. 



defining an elliptical annulus (distance bin), is: 
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Nmw and N are the expected number of MW and Scl RGB 
stars in a distance bin (iVx = Nmw+N). /mw is the velocity 
distribution of MW stars, which we assume does not change 
across the face of S cl, and is derived from the Besancon model 
( Robi net al.l 120031) selecting stars along the l.o.s. and with 
magnitudes and colors similar to the Scl RGB stars. We as- 
sume that the Scl velocity distribution is a Gaussian whose 
peak velocity v and dispersion a (the quantities we want to de- 
rive) are allowed to vary with projected radius. We derive the 
normalization factors, Nmw /Nt and N/Nt directly from the 
observed RGB surface density profile and relative foreground 
density. To estimate the fraction of MW interlopers in the MR 
and MP sub-samples we simply count how many stars with 
velocities < v sys — 3er (i.e. the non-membership region more 
populated by foreground stars) are classified as MR and as 
MP on the basis of their CaT derived [Fe/H] value. The like- 
lihood of observing a set of velocities Vi with i — 1, ...,N is 

L(vj, i>jv | v, a) = riiLi P( v i)- We maximize the likeli- 
hood function in each distance bin and find the corresponding 
best-fitting v(R) and a(R). The errors are determined from 
the intervals corresponding to 68.3% probability. 

The kinematics of the Scl MR and MP RGB stars are clearly 
different (Figure [3^,b): the l.o.s. velocity dispersion profile 
of MR stars declines from ^9 km s _1 in the center to ^2 
km s _1 at projected radius R = 0.5 deg, while MP stars are 
kinematically hotter and exhibit a constant or mildly declining 
velocity dispersion profile. 

4.2. Predicted Velocity Dispersion Profile 

The l.o.s. velocity dispersion predicted by the Jeans equa- 
tion for a spherical system in absence of net-streaming mo- 
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FIG. 3. — l.o.s. velocity dispersion profile (squares with errorbars), from 
rotation-subtracted GSR velocities, for the MR (a), MP (b) and all (c) RGB 
stars in Scl. The lines show the best-fitting pseudo-isothermal sphere (solid) 
and NFW model (dashed) in the hypothesis of j3 = /3om- Panel c) shows 
that the best-fitting pseudo-isothermal sphere with /3 = /3om (solid) and the 
NFW model with j3 =const (dashed) are statistically indistinguishable. 
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where R is the projected radius (on the sky), r is the 3D 
radius. The l.o.s. velocity dispersion depends on: the mass 
surface density £*(i£) and mass density p*(r) of the tracer, 
which in our case are the MR and the MP RGB stars; the 
tracer velocity anisotropy [3, defined as (3 = 1 — a 2 /a 2 , which 
we allow to be different for MR and MP stars; the radial ve- 
locity dispersion a r ,* for the specific component, which de- 
pends on the total mass distribution (for the general solution 
see Battaglia et al. l2005h . 

We consider two DM mass models: a pseudo-i sothermal 
sphere, typically core d, (see iBattaglia et al.l 120051). and an 
NFW profile, cusped dNavarro . Frenk & Whitd ll996). Since 
the contribution of the stars to the total mass of the sys- 
tem is negligible for reasonable stellar M/L ratios, we do 
not consider it further. As is unknown we explore two 
hypotheses: a velocity anisotropy constant with radius, and 
an Osipkov-M erritt (OM) velocity anisotropy (Osipko vl 19791 : 
Merrittl [l985l) . For the latter profile, the velocity anisotropy 
is j3 = r 2 j (r 2 + r 2 ) where r a is the anisotropy radius. 

4.3. Results from the Two-Components Mass Modeling 

We explore a range of core radii r c for the pseudo- 
isothermal sphere (r c = 0.001,0.05,0.1,0.5,1 kpc) and 
a range of concentrations c for the NFW profile (c = 
20,25,30,35). By fixing these, each mass model has two 
free parameters left: the anisotropy and the DM halo mass 
(enclosed within the last measured point for the isothermal 

8 We checked that the assumptions of sphericity and absence of streaming 
motions have a negligible effect on the results: the observed l.o.s. velocity 
dispersion profiles derived adopting circular distance bins and not subtract- 
ing rotation are consistent at the Icr level in each bin with the observed l.o.s. 
velocity dispersion profile derived adopting elliptical binning and by subtract- 
ing the observed rotation signal (see B07) 
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sphere, at 1.3 deg =1.8 kpc assuming a distance to Scl of 79 
kpc, see Mateo 1998; the virial mass for the NFW model). 
We compute a x 2 for the MR and MP components separately 
(Xmr an d XmP' respectively) by comparing the various mod- 
els to the data. The best-fitting model is obtained by mini- 
mizing the sum x 2 = Xmr + Xmp- We quote as errors in 
the individual parameters the projections of the A\ 2 = 3.53 
region (corresponding to the region of 68.3% joint probability 
for a three free parameters x 2 distribution). 

The models with constant velocity anisotropy yield poor fits 
for both a cored and a cusped profile. The best-fitting mod- 
els in the hypothesis of an OM velocity anisotropy are shown 
in Figure [5h,b and provide very good fits for both the MR 
and MP components, simultaneously. We find that a pseudo- 
isothermal sphere with r c =0.5 kpc, M(< 1.8 kpc) = 3.4 ± 
0.7 x 10 8 Mq, r a , M R = 0.21^5 kpc and r a , M p = 0.4±^ 
kpc gives an excellent description of the data (Xmin = 6.9, 
with Xmr = 0.6 and Xmp = 6-3). This gives an M/L ratio for 
Scl within 1.8 kpc of 158±33 (M/ L)^, an order of magnitude 
larger than previous estimates (lOueloz. Dubath & Pasquinil 
1995). Also an NFW model is statistically consistent with the 
data (c = 20, virial mass M v = 2.2±J;° x 10 9 M , r a;M R = 
0.2 ± 0.1 kpc and r a , M p = 0.8±^ kpc, X 2 min = 10 - 8 ' with 
Xmr = 4-2 and Xmp — 6.6), but tends to over-predict the 
central values of the MR velocity dispersion, and this ten- 
dency is accentuated for larger concentrations. This model 
gives a mass within the last measured point of ~ 2.4^' ^x 10 8 
M , which is consistent with the mass predicted by the best- 
fitting isothermal sphere model. 

Fig. fJJ; clearly shows the mass-anisotropy degeneracy 
present when modeling Scl as a single stellar component sys- 
tem: both a cored profile with OM velocity anisotropy and 
an NFW model with constant velocity anisotropy give a very 
good fit to the observed l.o.s. velocity dispersion profile, and 
they are indistinguishable. 



[Fe/H] measurements from VLT/FLAMES spectra of 470 Scl 
probable members. The large spatial coverage and statistics 
of this data-set, combined with the accurate velocity mea- 
surements, allowed us for the first time to detect in the stellar 
component of a dSph a statistically significant velocity gra- 
dient likely due to intrinsic rotation. This gradient is 7.6^2 '2 
km s _1 deg -1 along the Scl projected major axis. 

The presence of rotation in dSphs could support scenar- 
ios in which the progenitors of these objects were rotation- 
ally supported, disky sy stems, tidally stirre d by the interac- 
tion with the MW (e.g. iMayer et aill2001l) . The efficiency 
of the transformation from this kind of systems into dSphs, 
which are spheroidal in shape and have a pressure supported 
kinematics, is dependent on the eccentricity of the orbit and 
pericentric distance with respect to the host galaxy. It is likely 
that objects which have not been efficiently stirred could re- 
tain some of their initial kinematic characteristics, such as for 
example rotation. 

Under the hypothesis that Scl is not tidally disrupted we 
have determined its mass using the Jeans equations, which 
assume dynamical equilibrium. The combined fit of MR and 
MP stars allows us to to partly relieve the mass-anisotropy de- 
generacy present in modeling Scl as a single stellar population 
and to exclude the models with constant velocity anisotropy 
(Battaglia et al., in preparation). The best-fit is given by a 
cored profile of r c =0.5 kpc and mass M(< 1.8 kpc) = 
3.4 ± 0.7 x 10 8 M Q , assuming that the velocity anisotropy 
follows an O M model. This result is in a 1.5c agreement with 
the results of IStrigari et al.1 (120071) . who performed a careful 
comparison between different systems and datasets, measur- 
ing the mass of each galaxy within the same distance from 
the center. Our result makes Scl more massive than previ- 
ously thought, and this does not suppo rt the conjecture t hat 
all dSphs have a common mass scale dGilmore et al. 2007). 



5. DISCUSSION AND CONCLUSIONS 

We have explored the kinematic status and mass content of 
the Scl dSph using accurate line-of-sight velocities and CaT 



We thank L. Sales for kindly providing the orbits shown in 
Fig.l. 
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